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Abstract 

As is well known, the 0 — 0 component of the Schwarzschild space can be obtained by 
the requirement that the geodesic of slowly moving particles match the Newtonian 
equation. Given this result, we shall show here that the remaining components can 
be obtained by requiring that the inside of a Newtonian ball of dust matched at a 
free falling radius with the external space determines that space to be Schwarzschild, 
if no pathology exist. Also we are able to determine that the constant of integration 
that appears in the Newtonian Cosmology, coincides with the spatial curvature of the 
FLRW metric. 

1 Introduction 

General Relativity without Einstein’s equation 

It is interesting to see how far one can go in the derivation of GR results without using in fact the EE. 

Sommerfeld |1: tried to derive some aspects of General relativity starting from the Schwarschild metric 

ds 2 3 4 = - 7 ~ 2 dt 2 + 7 2 dr 2 + r 2 {d0 2 + si^Odcf 2 ) 

with 7 = 1 — AM-_ He then conceived that 


dto = 7 1 dt , 


(1) 


dr 0 = 7 dr , ( 2 ) 

He related them with time-dilatation and length-contraction. The time metric and space metric has 
been shown to be related with opposite power of 7 . But this construction has been shown to have many 
flaws [ 10 , P. Instead of pushing ideas along the lines of Sommerfeld, we investigate another direction, 
based on the fact that in cosmology, it is indeed possible to obtain the Friedmann equations for the 
expansion factor R(t). The cosmological solution may cut, thus representing a ’’ball ” of matter, outside 
this spherically symmetric ball we will assume there is a spherically symmetric static space time. Point to 
be noted that there is no use of Birkhoff’s theorem to provide the static assumption, here we have used 
the analogy of Newtonian gravity or electromagnetism where outside the matter or conducting sphere 
the field is static. As is well know, any static metric can be brought to the ’’standard” form, we will 
be labeling coordinates of the outside space using bars, so we define t and f for the time and radial 
coordinates outside. 


ds 2 = —B(r)dP + A(r)dr 2 + r 2 ( dO 2 + sin 2 0d(f 2 ) . (3) 

From the requirement that Newton’s law for gravity holds in some limit, we obtain an idea of what 
goo = —B should be: indeed, this metric component can be obtained by the requirement that the geodesic 
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of slowly moving particles match the Newtonian equation. For gff = A(f) is harder. 

In this paper we develop an approach which will go beyond the simple derivation of the 0 — 0 component of 
the Schwarzschild metric and give a derivation to obtain the g^f component under some assumptions that 
we will specify and avoid previous criticism [?]. The procedure consists of considering a free falling dust 
ball. As it is known that Newtonian cosmology [?] can be developed and a homogeneous dust system can 
be studied and by applying the Newtons laws.We can obtain that each dust particle expands or contracts 
according to a single expansion factor R, obeying 


87 tG 


k 

R 2 ’ 


( 4 ) 


where, k being an integration constant. This is possible for considering zero Cosmological constant. 
Points to be noted that this equation came out from the second order differential equation of R. 

So, the radius for a given particle is proportional to R{t), By using the continuity equation and the 
Euler motion equation, we can obtain that the density p ~ ^ and if we require -R(O) = 1 , i?( 0 ) = 0 , p = 
p( 0) at t = 0, we get 


k = 


87 tG 

— 


p{ 0) > 0 


and 


R(tf = k 



( 5 ) 


In the present work we will not require or derive or motivate a geometric interpretation of the constant 
of integration k. In contrast, a follow up section will see that, even with out using EE, k can be given 
the interpretation of spatial curvature in FRW space. Then matching to an external space can be used 
to determine the outside space by demanding continuity in some coordinate system from the inside to 
the outside. Of course, the face one can give a geometrical interpretation to the constant k without using 
EE is by itself interesting. Here we only will assume that the metric inside the ball is of the form 


ds 2 = —dt 2 + gijdx'dx^ 


( 6 ) 


and that the dust particles are co-moving that is dx l = 0 which are geodesics of the space time[?]. 
Furthermore, just from homogeneity and isotropy, we can determine that 


gijdx l dxi = R 2 {t) 


dr 2 

1 — nr 2 


+ r 2 {dd 2 + 



( 7 ) 


where the radial and time coordinates are r and t respectively. Notice that apriori we do not know if k 
and k are related. A very basic and elementary relation that we obtain when matching the space time 
© with ©-©, is obtained when considering the length of a circle of an equatorial circle (sinO = l)as </> 
runs from 0 to 27 t. The inside observer (in the dust ball), just below the surface of matching will measure 
a length 2 nrR, while The outside observer (in the static space), just above the surface of matching will 
measure a length 2 nf , since these two measurements refer to the same physical length, in order for the 
complete geometry to be well defined, we obtain 


r = rR(t). 


( 8 ) 


2 Newtonian Cosmology 

Following the Cosmological Principle, in a homogeneous and isotropic universe a spherical region of 
radius R can be identified around an arbitrary point, whose matter density p within is homogeneous. 
The surrounding matter cannot have any influence on its dynamics, as this would violate isotropy. By this 
general assumptions, the size of the sphere is arbitrary. The equation of motion for a test mass m located 
on the boundary of such a sphere shall be described in terms of a homogeneous positive parameter R(t), 
where the coordinate of each particle expands according to a(t) = constant ■ i?(t), where the constant 
depends on the particular particle. Along with the assumption of the motion of the particles as follows 
Bondi’s book [?]. 
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-$ = m-ct. 


( 9 ) 


P = P(t), 


( 10 ) 


P = P(t) = 0, 

here we are discussing the situation with P = 0. 

1 dR 


R dt 


= m ,R(t = o) = i, 


Now by applying Euler’s equation of motions , we have 

where the Hydrodynamic operator is + (u!V), 

now, taking divergence on both sides 


( 11 ) 


( 12 ) 


(13) 


3[f+ /’] = ?.* 


where, is the body force per unit mass. 
So by Poisson’s equation 


So, we have then 


7 .F = —AnGp 


3l f t+ f ] = - inGp 


(14) 


(15) 


(16) 


Putting the value of f(t) in the previous equation imply a similar equation for the universal expansion 
factor R(t) 


R=- 


G /4tt 
R 2 


R P = - 


47 tG 


-Rp. 


(17) 


Basically, this corresponds to Friedmann’s second equation without a cosmological constant A and 
zero pressure. As the linear dimensions scale by R(t), all co-moving volumes should scale by R(t) 3 , that 
is a 1/R 3 dependency for the density, which dilutes the matter as the sphere expands. 

For deriving a second equation, we first consider mass conservation within co-moving sphere, 


d ( 47t 
dt 


R P = 0, 


(18) 


where the internal mass M inside the sphere should be constant . By performing the derivative and 
simplifying one i?, the equation gets 


2 RRp + RRp + R~p = 0. 

The second term can be eliminated by m and after restoring derivatives the equation 

d(R 2 ) 8 t tG d{pR 2 ) 


dt 


dt 


(19) 


( 20 ) 


is obtained. Integration on both sides gives 


& = pR 2 _ k, 


( 21 ) 
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and rewriting the arbitrary integration constant k in a way to match the units k —> kc 2 yields finally 
an equation, which corresponds to Friedmann’s first equation in structure: 


R\ 

87tG ( i 

:\ 2 

u 

- 3 <• Hi 

?) 


( 22 ) 


It is interesting to point out the physical meaning of the constant k. Indeed, it is very much like a 
conserved energy of a mechanical system, and indeed the sign of k determines whether R will expand 
to infinity, or expand, achieve a maximum and then re collapse. Indeed, a simple analysis shows that 
for k > 0, we have a bound system, a bound system that has an associated negative energy, like wise 
k < 0 represents a case where the expansion factor reaches infinity and is the analogous of a mechanical 
system with positive energy. Finally, a case where R barely makes it to infinity is k = 0, and corresponds 
to a system with zero energy, analogous to the very similar problem of the critical escape trajectory 
that just makes in order to escape the earth. A separate issue is a possible geometrical interpretation 
of k 7 in particular a a possible relation between k and the k that appears in 0, this we will achieve, 
without assuming Einstein’s equations. The FLRW and the appearance of R(t) in 0 can be justified on 
purely symmetry and geometrical grounds, independent of any specific dynamics, General Relativity or 
otherwise. As we pointed out, although we will not invoke Einstein’s equations, we will make use of the 
much simpler geodesic equation. As we mentioned, it is a rather simple matter to show that ” co-moving 
observers”, of the metric ®, 0. that is those trajectories where the spatial coordinates in ©,0 are 
constants are indeed geodesics. Since, x l = const is a geodesic in FLRW spacetime as FLRW metric 
yields r£ = 0 [7] . Physical distances between these ” co-moving observers” scale as R , so the R in the 
Newtonian cosmology coincides with the R that appears in 0, the relation between k and n , this will 
require a bit more elaboration (without using EE). We can now describe cosmic large scale dynamics by 
a single expansion parameter R(t), as long the cosmic principle holds, and the predictions corresponds 
to GR. Note that from here on we will be using c = 1. 

3 Finding g 00 

The goo element of the metric is the first thing that students start form when learning GR, which a 
student might find hard to grasp. This is most likely the result of still not willing to be parted from 
Newtonian ideas, while accepting the ideas of special relativity. 

In this section, in order not to make the notation to heavy, we will for the moment denote the 
coordinates associated with the static space time simply t, r, 6 and <fi . In the next section, we go back to 
the bar notation when it refers to an external space and unbarred coordinates for the internal cosmological 
space time . In the next two subsections we will derive the gtt element of the metric from the geodesic 
equation in two different methods, which are really equivalent, but we want to present both to give the 
possible instructor some alternatives, or to present both. 

We allow our usage of the geodesic equation, because this subject does not involve a high level of 
complexity as compared to the development of the Einstein’s equations, and it can be derived from the 
from the principle of least action of the particle trajectory, the action being the proper time along the 
trajectory of the particle in a certain given metric. Starting from the interval 

ds 2 = g^dx^dx 1 ' = — dr 2 , (23) 

we integrate over dr and by using a monotonic variation parameter a on the particle trajectory 


dxv dx v 
9 ,iV da da 

from here applying the principle of least action we find the geodesic equation 

d 2 x 11 nU dx a dx^ 

+ °’ 



(24) 


(25) 
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with known as the Christoffel symbol (or connection) 

pM _ ( ^ 9 p 0 , 9 g ap __ dg a p \ 

a/3 2 y \ dx a dxP dxP J ' 


(26) 


Remember that is not a tensor, but the sum of the two terms in the geodesic equation is a vector. 

1st Method: Non - Relativistic limit 

e following method starts by a si 
part of the equation LHS, we have : 


The following method starts by a simple exchange of derivatives ^ = (£) 1 plugging in to the x l 


d 2 x l 

~dF 


-l 


d_ 

dr 


-l 


dx l 

dr 


dt\ 2 d. 2 x l 


dt \ 3 d 2 t dx l 


dr J dr 2 \ dr J dr 2 dr 


Now using the geodesic equation, the RHS can be re written as, 

” 2 rn ax ~ ax " , l at \ j 

a P dt dt dt 


dx a dx^ / dt 
a/3 dt dt + \dr 


3 ^0 da:“ dx^ dx l 


now equating both sides 


d 2 x i 


, dx a dx 13 0 dx a dx 13 dx 1 
dt 2 a/3 dt dt 0/3 dt dt dt 


(27) 


(28) 


(29) 


We continue this derivation by ordering in terms of v l = , by taking the following limit v 2 ~ hp and 

£=<?(!) 


dV 

dt 2 


= -n 


or* dxj r* dxj dxk 

dt jk dt dt 



dx 3 
1 dt 


t dxi dx k \ 

jk ~dT~dt J 


dx i 

dt 


(30) 


to derive the Newtonian limit we only consider the first term, since this is the term that dominates in 
the limit of very small velocities, 


= _ r i = _ 1 i P ( d 9 P t dg tp 
dt 2 u 2 y \ dx * dx * 


dgtt \ = 1 jj dg tt 
dx p ) 2 y dx? 


1 xij ®9tt 

2 dxi 


(31) 


In the above we have used that g lJ is in zero-th order approximation equal to 5 lJ , since the deviations 
from Minkowski space-time should be indeed very small. By demanding even for the small deviations to 
vanish at infinity, i.e. that asymptotic limit of the metric is exactly Minkowski space-time lim g tt = — 1 

r—>oo 

then by comparing with the Newtonian force law of universal gravity, 

d 2 x GM 

— T where </> =-, (32) 

dt- r 

we simply solve 

( 2 GM\ , 

9tt = — (1 ——) = — (1 + 2</>). (33) 


2nd Method: Metric perturbation 

The following method starts with taking two assumptions : 

1. The metric is static with a symmetric perturbation 

g P v — hp,v + hp.v j |^/xz/| 1 and 9pv ~ hp.u ? g^ ~ g k > (3d) 
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2. Non-relativistic limit 


, dx 1 „ dx 1 dt 

V = lit < ^ dT < d7 


We now return to the geodesic equation 


d 2 x^ ^„ dx* dx 4 „ da : 1 dad 

^ = “ r -»1717 = - rS * 17 ~ r « * 


- 1 - tt> 


writing the equations for the spatial part 

1 


F* — —-rT 

1 tt ~ n '/ 


%t\ 

Sad / 


= ^ 


/ dgtt 

V dx^ 


dV 

dt 2 


dgtt 

dx i 


(35) 


(36) 


(37) 


Now, by demanding the asymptotic limit lim g t t = —1, we once again comparing (21) to Newtonian 

r—>oo 

force law of universal gravity as in the previous sub-section 


9tt = - 1 - 


2GM ^ 


— — (1 + 2 4>). 


(38) 


As we mentioned before, the two methods are equivalent, but the first method gives us some idea of 
what terms we are dropping when doing our approximations, while the second method does not. 


4 Finding g ff 


The grr of the metric in the outside static region in our case, has been the ’’elusive component”, this 
component has not been calculated using matching to a Newtonian cosmology previously, here we will 
show that this is possible. In this section we will find gff using the assumption that we have a co-moving 
observer satisfying r = constant. We will also assume that inside and at the boundary of the dust ball 
the radius evolves as f = rR(t), where R(t) is determined by (l22l) . We also assume that for r > rR{t) 
the motion is on a radial geodesic and the metric is of the form 


ds 2 

dft 2 



2GM ^ 


dt 2 + A(f)df 2 + f 2 dn 2 


d6 2 + sin 2 9d<j ) 2 , 


(39) 

(40) 


we labeled the time as t because it may not be the same coordinate as the t in the dust ball and the gtt 
component was found in the previous section. Notice that we assume that the external metric is time 
independent. 

A radially falling geodesic, meaning that 9 = const and </> = const , is fully described by the conservation 
of energy that results from that the metric outside is assumed to be static. The geodesics are derived 
from the action 


S = 



dxv dx u . . 

do do 


The equation with respect to t is 


d / dL\ 
do V dt) 


0 


where 



This gives us 




2 GM'j 


dt 

dr' 


where 7 is constant and dr is the proper time. 

Notice that since the spatial coordinates in the space dx 1 = 0 we get 


dr = dt, 


(41) 


(42) 


(43) 


(44) 
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giving us 


dr 2 = 

(§) - - 


2 GM\ _ , 2 

1-— ) dt 2 - A{r)dr 2 


) 


r 

2 GM 


dr' 2 


Using previous equations, we obtain 


1 - 


2 GM 


7 -^-)(|) = 1 . 


(45) 

(46) 

(47) 


As we have seen, the consistency of the matching of the two spaces requires r = R(t)r , furthermore, 
we assume that even the boundary of the dust shell free falls according to a co-moving observer, which 
means that the FLRW coordinate r = constant and this allow then to solve for A(f), 


A(r) = - 1 - 


7 


We now use (3) and get 


1 


r 2 (^f) 2 ' 


(48) 


A(r) = - 1 - 


7 


1 


by simplifying this and expressing in terms of f we get 

Afr)- 1 7 2 -! + ^ 1 

l j .2 i _2 *(_!+£)• 


(49) 


(50) 


If we take the limit f —>• oo, we see that A(f) — (y 2 — l)/fcr 2 . Asymptotic flatness would require 
A(f) —> 1. 

The metric component A(r) is free of singularities (real singularities, not coordinate singularities) and 
preserves its signature (one time and three spaces) and is asymptotically flat only if 


y 2 — 1 = —kr 2 and 


2 GM r 2 k kr 3 


R 


k = 


2 GM 


r U r r° 

As mentioned earlier, r is the comoving coordinate and hence constant, whereas R is variable. 


(51) 


Notice that above the condition. (1511 k = 2G t I , when combined with the value of k, as given by ©, 


k = sikp(O) yields us the value of M as 


M = —TTp 0 r 3 . 


(52) 


Surprisingly, it looks like the Mass = density~x.Vol.ofE 3 ball. But, we got this relation in consequences 
of our previous derivation and the flat space volume of the ball has not been used anywhere for deriving 
this. 


Finally, all of this gives us 


A{r) = 


1 


1 - 


2 GM ’ 


(53) 


reproducing the Schwarzschild spacetime. 
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5 Finding the geometric interpretation of the Newtonian Cos¬ 
mology 

The main aim of this section is to show that for the geometric parameter defining the FRW space 
coincides with the Newtonian Energy in k found from the interpretation of the Newtonian cosmology and 
to demonstrate that k = k. In an embedding four dimensional Euclidean Space with metric 

dl = dx T dy T dz T dw , 

we define a 3-sphere if the sets of points ( x, y , z , w) satisfy the constraint 

x 2 + y 2 + z 2 + w 2 = 1 /k with n > 0, 

then it is said to form a three - sphere S 3 . That is, we will take our cosmological solution to have positive 
spatial curvature k > 0 therefore topologically spatial slices have the topology of S 3 . Solving for w from 
the above constraint and inserting into the expression for dl 2 defining r 2 = x 2 + y 2 + z 2 , x = rsinOcosrj) 
and y = rsin9sincj),z = rcosQ , we obtain the metric of the 3-sphere, 

ds l = [ , dr 2 +r 2 dn 2 ], 

1 — nr 2 

with 

dfl 2 = {dd 2 + sin 2 0d(j) 2 ), 

to obtain the physical space time, that is the line element in FLRW Space, we add time and a scale factor 
that multiplies the 3-sphere, we obtain, 


ds 2 = —dt 2 + R 2 (t).ds\. 

So our infinitesimal line element in FRW space becomes, 


ds 2 = —dt 2 + R{t) 2 ( — — + r 2 dtt 2 ') . 

\ 1 — ) 


(54) 


We assume again that that we have a co-moving observer which satisfies r = const. Independently of 
that, in the FLRW space we can use everywhere (not just at the boundary) the barred radius f = R(t)r , 
which means r = 

This yields 


RW 


dr Rr 
dr = -- Trdt. 

R(t) R 2 


Putting this in (35), we’ve 


(55) 


ds 2 = -1 + 


R 2 f 


(1 - nr 2 )R 2 


dt 2 


di' 2 


2 Rf 


1 — nr 2 (1 — nr 2 )R 


dtdr. 


Now, we make a transformation t = t(t,r), so infinitesimal change in time 

, dt dt - 
dt = —dr + -- dt , 
dr dt 


squaring, 


dt 


dr 


dt 


dt 


dt 2 = — dr + ( —- ) dt + 2 — -= dtdr. 


dt\ (dt 


dr \dt 


(56) 


(57) 


(58) 


Under these circumstances, the infinitesimal line element becomes 



(59) 


ds 2 = 1—1 


R 2 r 


(1 - Kr 2 )R 2 


dt 


dr 


dt 


dt 


dr + -= dt + 2 — dtdr 


dt\ (dt 


+ 


dr 2 


2 Rr 


dr j \dt 


, , dt dt ' 
dr — dr + —zdt 


1 — nr 2 (1 — nr 2 ) R \dr dt 


Now, we have to eliminate the cross terms, so 


2 Rr dt - ( R 2 r 

-dr-^dt+ -1 + 


This gives 


(1 — nr 2 )R dt 


(1 — nr 2 )R 2 ) \dr J \dt 


dt\ (dt 


dtdr = 0 . 


rR 


(1 — nr 2 )R ^ 


_i_ 

1 + (1 —Kr 2 )i? 2 


rR 


R 2 r 


-(1 - nr 2 )R H— 
fRR 

-(1 - nr 2 )R 2 + R 2 f 2 ' 


Now lets linger on g f f with this new (44) value 


1 


Qrr — 


1 — nr 2 

1 

1 — Kr 2 

1 

1 — Kr 2 


+ -1 


+ -1 


+ -1 


2 Rr 


R 2 r 


2=2 


i dt 2 2 Rr dt 

(1 — nr 2 )R 2 ) dr (1 — nr 2 )R dr 


R 2 r 2 


dt 


2 Rr dt 


(1 — nr 2 )R 2 j \dr J (1 — nr 2 )Rdr 


R 2 r 2 


R 2 r 2 


(1 - nr 2 )R 2 J _ Kr 2) + 
Rf 


R 2 f 


T )W 


(l-^ 2 )«i!(l-^)(-l + TT Jg^) 


1 — nr 2 


1 + 


R 2 f 2 


2 R 2 f 2 


-R 2 + tiR 2 r 2 + R 2 f 2 -R 2 + nR 2 r 2 + R 2 f 2 / 
R 2 ( 1 — Kr 2 ) 


(1 — Kr 2 ) (^—R 2 + kR 2 t 2 + R 2 r 2 ^j 

1 


1 - Kr 2 - 


For r = R(t)r it yields 


(60) 

(61) 

(62) 

(63) 

(64) 


Qrr — 


1 

1 — r 2 (ft + R 2 ) 


In the last section we derived also what this metric component should be, 

1 

9rr ~ ^ 2 GM ' 

f 


(65) 


( 66 ) 
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Since from Newtonian cosmology , we know that 


(67) 


Therefore, 


k + R 2 


8n 

3R 


G Po 


8irr 3 p 0 
3r 2 (rR)' 


(.k + R 2 )r 2 


8nr 3 p 0 
3 Rr 


2 GM 
f 


( 68 ) 


recalling that M = ^irpor 3 and f = rR giving us the relation k = n . 

A final consistency check is obtained now that we have derived that the internal space is a Robertson 
Walker space, satisfying the standard Friedmann equations, with k being interpreted as the spatial 
curvature, as we have just shown. Under this condition, we know the matching of this internal cosmology 
to Schwarzschild is consistent, as the analysis of the Oppenheimer- Snyder collapse model shows[jQ. Here 
we have gone about this problem in the opposite way, showing that the matching of these two spaces 
imposes severe constraints that allows us to derive Schwarzschild space in the outside and determine that 
the Newtonian constant of integration k has to be the spatial curvature of the internal FLRW internal 
space, all of this without using Einstein’s equations. 


6 Conclusions 

In this paper it has been found that matching a dust ball, whose dynamics is governed by the Newtonian 
cosmology equations, containing a constant of integration fc, to an external static space-time, where 

ds 2 = —B(r)dP + A(r)df 2 + r 2 dQ 2 ( 69 ) 

where B = (l — 2G ~ M ), forces A to have a very special form. Assuming either asymptotic flatness, or 
absence of signature change, we uniquely obtain 

^( r ) = ^ 2GM~ y ’ (70) 

Finally the same matching to the internal FLRW space 

ds 2 = — dt 2 + R 2 (t) ( ——— -r + r 2 dQ 2 '\ ( 71 ) 

\ 1 — k r A ) 

forces the geometrical parameter k that appears in FLRW to coincide with the constant of integration k, 
used in the Newtonian Cosmology. 

In this essay, it has been found that matching a dust ball, whose dynamics is governed by the Newto¬ 
nian cosmology equations, containing a constant of integration k coincides with the geometrical parameter 
k that appears in FLRW used in the Newtonian Cosmology. These results are of interest at least in two 
respects, one from the point of view of its pedagogical value of teaching General Relativity without in 
fact using Einstein’s equation and second, the fact that some results attributed to General Relativity 
can be obtained without using General Relativity indicates that these results are more general than the 
particular dynamics specified by General Relativity. Although, some generalizations are possible, like 
the possibility of introducing a cosmological constant in Newtonian Cosmology, as discussed in Bondis 
book [ 5 ] , we could in this way that by matching this cosmology to an exterior stationary space we obtain 
Schwarzschild- deSitter space, but EE do not expect that an approach of this type will be able to give 
all results of GR, for example, certainly not in the case of gravitational waves or the Kerr solutions.For 
the case of inhomogeneous dust ball distribution, the CM lies on arbitrary origin and hence the linear 
and angular momentum remains zero. So also has been shown that the effective force between the con¬ 
stituent particles are zero for no perturbation [S] .If we ponder minutely, the paper [5] goes very far just 
using Newtonian physics, but is missing a space-time interpretation of the solutions of Newtons laws for 
gravitationally interacting particles. We might try to find them properly even in presence of pressure in 
our next venture. 
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Finally an interesting question arises that whether our derivation holds only in the weak field approx¬ 
imation or not. Notice that indeed in parts of our arguments we have used the weak field approximation, 
like when we derived the 0 — 0 component of the metric, then on this basis, we derived the r — r component 
by a process of matching to the internal collapsing ball of dust. But if we take the point of view that 
we trust the metric of the collapsing dust beyond the weak field approximation, the situation will be 
different, in this case our derivation could have validity beyond the weak held approximation, this is a 
question to be studied. One should also point out that in GR, as far as the post-Newtonian approximation 
is concerned, the corrections to the the 0 — 0 component of the metric appear at the same other as the 
first corrections to the r — r components, so it is in a sense puzzling that the r — r component has been 
more elusive to find by a simple derivation, as compared to the case of the 0 — 0 component. 
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